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Abstract 

The Lamperti-Kiu transformation for real-valued self-similar Markov processes 
(rssMp) states that, associated to each rssMp via a space-time transformation, there 
is a Markov additive process (MAP). In the case that the rssMp is taken to be an 
a-stable process with a € (0,2), [16] and [23] have computed explicitly the charac¬ 
teristics of the matrix exponent of the semi-group of the embedded MAP, which we 
henceforth refer to as the Lamperti-stable MAP. Specifically, the matrix exponent of 
the Lamperti-stable MAP’s transition semi-group can be written in a compact form 
using only gamma functions. 

Just as with Levy processes, there exists a factorisation of the (matrix) exponents of 
MAPs, with each of the two factors uniquely characterising the ascending and descend¬ 
ing ladder processes, which themselves are again MAPs. To the author’s knowledge, 
not a single example of such a factorisation currently exists in the literature. In this 
article we provide a completely explicit Wiener Hopf factorisation for the Lamperti- 
stable MAP. 

The main value and novelty of exploring the matrix Wiener Hopf factorisation of 
the underlying MAP comes about through style of the computational approach. Un¬ 
derstanding the fluctuation theory of the underlying MAP offers new insight into dif¬ 
ferent ways of analysing stable processes. Indeed, we obtain new space-time invariance 
properties of stable processes, as well as demonstrating examples how new fluctuation 
identities for stable processes can be developed as a consequence of the reasoning in 
deriving the matrix Wiener-Hopf factors. The methodology in this paper has already 
lead to new applications in the forthcoming work of m and [28] . 

Key words: Self-similar Markov process, Lamperti-Kiu transform, Markov additive 
factorisation, matrix Wiener-Hopf factorisation. 
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1 Introduction 

Let X := (A5)t>o be a one-dimensional Levy process, starting from zero, with law P. The 
Levy-Khintchine formula states that, for all 6 6 M, the characteristic exponent \l/(#) : = 
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logE(e lPXl ) satisfies 


\P(0) = iad + ^rcr 2 d 2 + f (1 — e l8x + i9xl(\ x \<i))U(dx), 9 G M, (1) 

^ Jr 

where a G M, a > 0 and II is a measure (the Levy measure ) concentrated on M\{0} such that 
J R (1 A x 2 )II(dx) < oo. When analytical extension is possible, we refer to ijj(z) := — \P(— iz) 
as the Laplace exponent. 

The process (A", IP) is said to be a strictly a-stable process (henceforth just written ‘stable 
process’) if it is an unkilled Levy process which also satisfies the scaling property, under P, 
for every c > 0, the process ( cX tc - a ) t > 0 has the same law as X. It is known that a necessarily 
belongs to (0,2], and the case a = 2 corresponds to Brownian motion, which we exclude. 
The Levy-Khint chine representation of such a process is as follows: a = 0, II is absolutely 
continuous with density given by 

vr(x) := c+x~ {a+1) 1 ( X>0 ) + c_|x|” (q+1) 1 (3; <o), x E R, 

where c+, > 0, and a = (c+ — c_)/(a — 1). 

The process X has the characteristic exponent 

V(G) = c\9\ a (l - i/3tan^sgn(0)), 9 e M, (2) 

where f3 = (c + — c_)/(c+ + c_) and c = — (c+ + c_)T(— a) cos(7ra/2). Self-similarity dictates 
that we must necessarily have /3 — 0 when a — 1, which is to say that the process is 
symmetric. For more details, see [33j Theorems 14.10 and 14.15]. 

For consistency with the literature that we shall appeal to in this article, we shall always 
parametrise our cc-stable process such that 

sin(7rap) , lN sin(7ra;p) 

c + = T(a + 1)—-—— and c_ = r(a + 1)—-——, 

7r 7r 

where p = P(A t > 0) is the positivity parameter, and p = 1 — p. In that case, the constant 
c simplifies to just c = cos(7ra(p — 1/2)). Moreover, we may also identify the exponent as 
taking the form 

*(0) = |6'|“(e 7ria( 5 _p) l (0 >o) + e-™^- p )l (0<o) ), 9 e M. (3) 

With this normalisation, we take the point of view that the class of stable processes is 
parametrised by a and p\ the reader will note that all the quantities above can be written 
in terms of these parameters. We shall restrict ourselves a little further within this class by 
excluding the possibility of having only one-sided jumps. In particular, this rules out the 
possibility that A is a subordinator or the negative of a subordinator, which occurs when 
a 6 (0,1) and either p = 1 or 0. In the case of a subordinator, for future reference, we note 
that the characteristic exponent takes the form T(6 I ) = (—i 9) a , 9 e M, which is the analytic 
extension of the Bernstein function A i —> A“, A > 0. 

A fascinating theoretical feature of all characteristic exponents of Levy processes is that 
they can always be written in terms of the so-called Wiener-Hopf factors. That is to say, for 
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a given characteristic exponent of a Levy process, T, there exist unique Bernstein functions, 
k and k such that, up to a multiplicative constant, 


^( 0 ) = k(i9)n(—i9), 9 G 


(4) 


As Bernstein functions, k and k can be seen as the exponents of (killed) subordinators. The 
probabilistic significance of these subordinators, known as the ascending and descending 
ladder height processes respectively, is that their range corresponds precisely to the range of 
the running maximum of X and of —X respectively. In this sense, they play an important 
role in understanding the path fluctuations of the underlying Levy processes. In particular, 
a rich history of literature has shown their fundamental significance in the understanding of 
a variety first passage problems; see for example their extensive use the the development of 
fluctuation theory of Levy processes in the texts |8] and [26] as well as m. 

In the case of stable processes, the Wiener-Hopf factorisation takes a relatively straight¬ 
forward form. Indeed, it is straightforward to argue that the ascending and descending ladder 
processes must necessarily be stable subordinators. One is therefore forced to take (up to a 
multiplicative constants) /c(A) = A ai , A > 0, and k( A) = A Q2 , A > 0, for some ai ,«2 € (0,1). 
Comparing ([4]) with (3), we must choose the parameters aq and aq such that, for example, 
when z > 0, 


z e 


7ria(|-p) _ 


= z 


« 1 „- 


2 Trial ^ 


(5) 


Matching radial and angular parts, we find that 


oq T 0-2 — cr, 

aq — «2 = — a(l — 2 p), 


( 6 ) 


which gives us aq = ap and «2 = otp. As we have assumed that X does not have monotone 
paths, it is necessarily the case that 0 < ap < 1 and 0 < ap < 1. Note also that when 
ap = 1, the ascending ladder height process is a pure linear drift. In that case, the range 
of the maximum process X is [0,oo). This can only happen when X is spectrally negative 
which has been ruled out by assumption in the introduction. Similarly the case that ap — 1 
corresponds to spectral positivity which has also been ruled out by assumption. In conclusion, 

k(A) = X ap and k( A) = X ap , A > 0 


where 0 < ap, ap < 1. 

This discussion also helps us justify that, taking account of all the special cases of stable 
processes that we have chosen to exclude, the set of admissible parameters we are left to 
work with is 


{(a, p) : a G (0, 2), p G (1 — 1/a, 1/a) and p=l/2ifa = l}. 

In this article, we expose a second Wiener-Hopf factorisation which is ‘deeply’ embedded 
within the stable processes through its so-called Lamperti-Kiu representation. Like the 
factorisation (|5|), the ‘deep factorisation’ we will present has value in that it informs us 
about the fluctuations of the stable process. As we shall see, this gives us access to an array 
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of new results for stable processes, as well as a methodology for obtaining even more than 
those presented in this paper. 

The Lampert-Kiu representation is a pathwise decomposition that holds more generally 
for any real-valued self-similar Markov processes (rssMp) and shows that any such process 
can be written as a space-time changed Markov additive process (MAP). In a similar spirit to 
©. the semi-group of a MAP can be characterised via an exponent, albeit that it now takes 
the form of a complex-valued matrix function. Moreover, just as with Levy processes, there 
exists a factorisation of the aforesaid matrix exponent. Not a single concrete example of such 
a factorisation currently exists in the literature for such MAPs to the author’s knowledge. 
Our main objective here is to provide a completely explicit Wiener-Hopf factorisation for 
the MAP that underlies the stable process via the Lamerti-Kiu transform: the so-called 
Lamperti-stable MAP. 

Our approach appeals to three main techniques. First, information about the the respec¬ 
tive individual entries in the matrix factors can be gleaned using asymptotic Markov additive 
renewal theory in the setting of excursion theory for MAPs. Second, the quantities that are 
identified in that way can be related to complex first passage problems for stable processes 
Third, the aforesaid first entry problems can be simplified by appealing to a version of the 
Riesz-Bogdan-Zak transform, which relates the mapping of the path of a stable process via 
a Kelvin transform, together with an endogenous time change, to a Doob /i-transform of the 
stable process. For a G (1, 2), this Doob h-transform corresponds to conditioning the stable 
process to avoid the origin, as explored in For a G (0,1) it corresponds to conditioning 
the stable process to being absorbed at the origin. Finally, for a = 1 there is, in effect, no 
h-transform as h = 1. 

The remainder of this article is structured as follows. In the next section we explain 
the nature of the Lamperti-Kiu representation for pssMps, due to [TB] , and its relation to 
MAPs. In particular we give the example of the of the Lamperti-stable MAP, also due to 
[16]. With this in hand, we are able to state our ‘deep’ Wiener-Hopf factorisation result 
for the stable process. In section [3] we discuss the Riesz-Bogdan-Zak transform. In section 
0 we compute the first matrix factor of the Wiener-Hopf factorisation for the Lampert- 
stable MAP, which is the analogue of the contribution from the exponent of the ascending 
ladder height process in the Wiener-Hopf factorisation for Levy processes. In section [7] we 
compute the second Wiener-Hopf matrix factor, which is the analogue of the contribution 
from the ascending ladder height process of the dual in the Wiener-Hopf factorisation for 
Levy processes. Finally, in Section [SJ we outline how some Cramer-type asymptotics and 
explicit identities for the Lamperti-stable MAP and stable processes can be obtained from 
the methods that underly the deep factorisation. 

2 MAPs and the Lamperti-Kiu transform 

This section is laid out as follows. We devote the first two subsections to a discussion 
of Markov additive processes and real self-similar Markov processes via the Lamperti-Kiu 
representation. Finally, in the last subsection, give our main result, the deep Wiener-Hopf 
factorisation of the stable process. 
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2.1 Markov additive processes 

Let E be a finite state space and (Gt)t> o a standard filtration. A cadlag process (£, J) in 
Kx£ with law P is called a Markov additive process (MAP) with respect to (Gt)t> o if 
o is a continuous-time Markov chain in E , and the following property is satisfied, for 
any i E E, s,t > 0: 


given { J(t) = i}, the pair (£(t + s) — £(f), J(t + s)) is independent of Gt, 
and has the same distribution as (£(s) — £(0), J(s)) given {J(0) = i}. (7) 


Aspects of the theory of Markov additive processes are covered in a number of texts, 
among them [6j and [7||. We will mainly use the notation of [T9] , where it was principally 
assumed that £ is spectrally negative; the results which we quote are valid without this 
hypothesis, however. 

Let us introduce some notation. For i£l, write P Xji = P(-1 £(0) — x, J( 0) = i). If /x is 
a probability distribution on E, we write P I)AJ = JA gE HiP x ,i- We adopt a similar convention 
for expectations. 

ft is well-known that a Markov additive process (£, J) also satisfies (7) with t replaced by a 
stopping time, albeit on the event that the stopping time is finite. The following proposition 
gives a characterisation of MAPs in terms of a mixture of Levy processes, a Markov chain 
and a family of additional jump distributions; see [Tj §XI.2a], |T9j Proposition 2.5] as well as 
more classical literature such as mm®- 


Proposition 1. The pair (£, J ) is a Markov additive process if and only if, for each i,j E E, 
there exist a sequence of iid Levy processes (£f) n >o an d a sequence of iid random variables 
{Ufj) n > o, independent of the chain J, such that if cr 0 = 0 and (cr n ) n >i are the jump times of 
J, the process £ has the representation 


£(t) = l (n> o)(f(<Tn-) + f "(*„)(*“ O > * ^ [<T n ,CT n+1 ), 71 > 0. 

For each i e E, it will be convenient to define, on the same probability space, £j as a Levy 
process whose distribution is the common law of the £” processes in the above representation; 
and similarly, for each i, j G E, define to be a random variable having the common law 
of the Uf - variables. 

Henceforth, we confine ourselves to irreducible (and hence ergodic) Markov chains J. 
Let the state space E be the finite set {1,. . . , A/”}, for some N G N. Denote the transition 
rate matrix of the chain J by Q = ( Qi,j)i,jeE■ F° r each i G E, the Laplace exponent of the 
Levy process £,; will be written For each pair of i,j G E , define the Laplace transform 
Gij(z) = E(e zUi ' j ) of the jump distribution U t ^, where this exists. Write G(z) for the N x N 
matrix whose (i,j) th element is G l j (z). We will adopt the convention that U i} j = 0 if q ltJ = 0, 
i ^ j, and also set Ua = 0 for each i G E. 

The multidimensional analogue of the Laplace exponent of a Levy process is provided by 
the matrix-valued function 


F(z) = diag(V»i(«), • ■ ■, Vtv(a)) + Q ° G(z), 


( 8 ) 
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for all z G C where the elements on the right are defined, where o indicates elementwise 
multiplication, also called Hadamard multiplication. It is then known that 

Eo,j(e'«‘>; J(t) = j) = (e F W‘). i,jeE,t> 0, 

for all ^gC where one side of the equality is defined. For this reason, F is called the matrix 
exponent of the MAP (£, J). 

Just as is the case with Levy processes, the exponents of MAPs are also known to have 
a Wiener-Hopf factorisation. However, this time, the two factors correspond to the matrix 
exponent of the ascending (resp. descending) ladder processes. These are themselves MAPs 
with trajectories which agree with the range and state of the modulating chain at times of 
new maxima (resp. minima). In order to explain the nature of the Wiener-Hopf factorisation 
for MAPs, we need to introduce a little more notation. 

Associated to the running maximum process (sup s<t £(s))t>o is a Markov additive subor- 
dinator. That is, a MAP, say ( H + {t ), J + {t )) t > 0 , with the property that H + is non-decreasing 
with the same range as the running maximum. Moreover, its exponent can be identified by 
—k(—z), where 

«(A) = diag($ 1 (A), • • • , <MA)) - A o K( A), A > 0, (9) 

is a matrix analogue of a Bernstein function. Here, for i = 1, • • • , N, <f>j are Bernstein 
functions (exponents of subordinators), A = (A i,j)i,jeE is the intensity matrix of J + and 
= E[e- w S], where f/T > 0 are the additional discontinuities added to the path of 
£ each time the chain J + switches from i to j, and C/L := 0, i G E. 

We also need to talk about the same quantity but for the dual of (£, J). Whilst the 
dual of a Levy process is equal in law to nothing more than its negative, the situation for 
MAPs is a little more involved. First note that, thanks to irreducibility, the Markov chain J 
necessarily has a stationary distribution. We denote it by the vector 7V — (tt 1 , ■ ■ ■ The 

dual process that is the MAP with probabilities P Xj j, x G M, i e E, whose matrix exponent, 
when it is defined, is given by, 

E„, J(f) = j] = (e FW ‘). j, t, j 6 E, 

where 

F(z) := diag (i/j^-z), ...,ip\ E \(-z)) + Q o G(-z) T 
and Q is the intensity matrix of the modulating Markov chain on E with entries given by 

7 Tj 

Qij = — Qj,i, i,jeE. 

Note that the latter can also be written Q = A" 1 Q T where A,r = diag(7r), the matrix 
with diagonal entries given by 7r and zeros everywhere else. Hence, when it exists, 

F(z) = a; 1 f(- 2 ) t a„ (io) 

showing that 

vqEo,* [e^W, J(t) = j] = 7t,E 0j J(t) = i ]. (11) 

At the level of processes, one can understand (fill) as saying the following. 
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Lemma 1. The time-reversed, process {(£((f — s)—) — £(i), J((t — s )—)) : s < t} under Po.^ 
is equal in law to {(£(s), J(s)) : s < t} under P 0i7r . 

We are now ready to state the Wiener-Hopf factorisation for MAPs. Whilst some results 
in this direction exist in classical literature, see for example Chapter XI of [5] or Theorem 
3.28 of [21]. Several other references can be cited in this respect, for example [23] and [4], 
None of them are in an appropriate form for our purposes. We have lifted the following 
result from the Appendix of the recent article [fB]. 

Theorem 1. For 6 G R, up to an multiplicative constant, 

-F(i6) = A^k^A^-W), 


where k plays the role of k, but for the dual MAP to (£, J). 


Note that this Theorem is consistent with the Wiener-Hopf factorisation for Levy pro¬ 
cesses 


(4) as, in that setting, the dual process is its negative. 


2.2 Real self-similar Markov processes 

The structure of real self-similar Markov processes has been investigated by on in the 
symmetric case, and [T6] in general. Here, we give an interpretation of these authors’ results 
in terms of a two-state Markov additive process. We begin with some relevant definitions 
and introductory results. 

A real self-similar Markov process (rssMp) with self-similarity index a > 0 is a standard 
(in the sense of na ) Markov process X = (X t ) t>0 with probability laws (P x )a:eR which 
satisfies the scaling property that for all i6l \ {0} and c > 0, 

the law of (cX tc -<x) t > 0 under is P cx . 

In [IB] the authors confine their attention to processes in ‘class C.4’. An rssMp X is in 
this class if, for all x ^ 0, P x (3t > 0 : X t X t _ < 0) = 1; that is, with probability one, the 
process X changes sign infinitely often. Define 

r {0} = inf(t > 0 : X t = 0}, 

the time to absorption at the origin. 

Such a process may be identified with a MAP via a deformation of space and time which 
we call the Lamperti-Kiu representation of X. The following result is a simple corollary of 
[161 Theorem 6]. 

Proposition 2. Let X be an rssMp in class C-4 and fix x ^ 0. Define the symbol 

1, y > o, 

2, y< 0. 

Then there exists a time-change a, adapted to the filtration of X, such that, under the law 
P x , the process 

m,j(t)) = ( logix^i,^]), 

7 


t > 0, 



is a MAP with state space E = {1,2} under the law Pio g p|,[x]- Furthermore, the process X 
under P x has the representation 


X t = exp(f(<p(t)) + m(J(ip(t)) + 1)), 0 < t < r {0} , 

where <p is the inverse of the time-change cr, and may be given by 

ip(t) — inf is > 0 : [ exp(a£(w)) dw > t\, t < t®, 


( 12 ) 


such that (£, J) has law Piogxjx]- 


2.3 The Lamperti-Stable MAP and deep factorisation 

Now let us return to the case that X is a stable process as described in the introduction, which 
is also an rssMp. In pT6l §4.1], the authors calculate the characteristics of the Lamperti-Kiu 
representation for X until absorption at the origin; that is, they compute the characteristics 
of the processes £*, the jump distributions Uij and rates the q it j, for i,j G {1, 2}. Using this 
information, and the representation |(8)[ it was shown in m that the MAP (£, J) has matrix 
exponent 


F(z) 


r (o? — z)r(i + *) 

Y(ap — z)r(l — ap + z) 

r(a — z)r(i + z) 
r(ap)r(l — ap) 


r(a — z)r(i + z ) 
r(ap)r(l — ap) 

r(a — z)T(l + z) 
r(ap — z)T(l — ap + z) 


(13) 


for Re(z) G (—1, a). In the spirit of [2B, Chapter 13.4] we refer to this process as a Lamperti- 
stable MAP. 

We should also note that the diagonal terms have entries which are characteristic expo¬ 
nents which belong to the class of so-called hypergeometric Levy processes. Moreover, up to a 
multiplicative constant, the off-diagonal terms can be shown to be the Laplace transforms of 
distributions, which possess a density with respect to Lebesgue measure that can be written 
in terms of the classical hypergeometric 2 X\ function (see Chapter 13 of [2B]). The matrix 
exponent (fl3]i could, in theory, be shown to belong to a bigger family of MAPs which are 
in some sense a natural generalisation of the class of hypergeometric Levy processes (cf. 
|23]). Indeed, we shall see other MAPs in the forthcoming analysis which are different to the 
Lamperti-stable MAP but clearly are close relatives with a common analytic structure. We 
shall explore this remark in more detail in future work however. 

Our main result, below, gives the explicit factorisation of (|T3]) as predicted by Theorem 
[U To our knowledge this is the first time that an example of the Wiener-Hopf factorisation 
for a MAP has has been detailed explicitly. 


We first need to introduce some notation. Of use will be the family of Bernstein functions 


^q+itP+ji^) 


(1 - e~ Xx ) 


((q + i) V(p + j)-l) ax 
(1 - e- x ) 9+i (l + e~ x )P + j 
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o 


A > 0, (14) 











where q,p G {ap, ap} and i,j G {0,1} such that q + p = a and i + j = 1. Note that it is easy 
to verify that the above expression is indeed a Bernstein function as the associated Levy 
density behaves like either x _ " p_1 or x _ “ p_1 as x j. 0 and like e~ ax as x f oo. Accordingly, 
it is also straightforward to verify that the mean value /c^ +ip+J -(0+) is finite. 

Theorem 2. When a G (0,1] , We have the following two components to the factorisation 
in Theorem[ 7J 


(i) Up to a multiplicative constant, the ascending ladder MAP exponent is given by 






Q!/9+l,Q'P 


Qp(^)' 


s'm(irap) 
sin(7 rap) 


exp, ap+1 


(o+) 


sin(Trqp) « n p,ap+l(N 
sin(7rap) A 


sin(7rap) K ap , n p +1 (A) 
sm(nap) A 


l^exo- 


cxp-\-l,cxp 


(A) 


sin(7rap) , 
sin(Trap) K <*P’ a P+ 1 


( 0 +) 


for A > 0. 

(ii) Up to a multiplicative constant, the dual ascending ladder MAP exponent is given by 


k(X) = 


sin(7rc*!p) / 

l.a/.fA + 1 - «) + sinW) %,aJ+l 
K a/5,ap+l(^ + 1 — a) 

A + 1 — Q! 


( 0 +) 


K ap,Qp+l(^ “t - 1 01 ) 

A + 1 — a 

sin(7ra:/3) / 

^'Q'P+ljQ'pv^ H - 1 “f - ; 7 r U’n 


sin(7r ap) a ^ +1 


( 0 +) 


for A > 0. 


The next theorem deals with the case that a 
another family of Bernstein functions. Define 


G (1,2). For this we need to introduce 


< i>q+i,p+j{^) ^ (1 


—Xu\ 


((q + i ) v (p + j) - 1 ) 


(« - 1 ) 


(i 


l Y +i ( 1 + e~ u )p+o 2(1 - e _u )' ? (l + e~ u )P 


Mu, 


for A > 0, q,p G {ap, ap} and i,j G {0,1} such that q + p = a and i + j — 1. Note, again, 
that the density in curly brackets can easily be verified to be positive in all cases and is a 
Bernstein function since, as before, the associated Levy density behaves like either x _ap_1 
or x _ “ p_1 as x 0 and like e~ x as x f oo. Once again, it is also subsequently straightforward 
to verify that the mean value /c' +ip+J -(0+) is finite. 

Theorem 3. When a G (1,2), we have the following two components to the factorisation in 
Theorem QJ 


(i) Up to a multiplicative constant, the ascending ladder MAP exponent is given by 


k(A) = 


sm(irap)<f> ap+1}a p(\ + a - 1) + sin(wap)t/>' a ^ ap+1 (0+) 

Pap,ap+lU + 0 — 1) 


— sin(7TQ!/5) 


0o:p,Q:p+ l(A + Q! 1) 


— sin(7rap)- 


A + a — 1 


A ex, — 1 

sm(nap)<t> a p +1 , ap (\ + a - 1) + sin(irop)^ pa/3+1 (0+) 


for A > 0. 
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(ii) Up to a multiplicative constant, the dual ascending ladder MAP exponent is given by 


k(X) = 


sin(irap)<j> a p +l!Cip (X) + sin^ap)^ ^ +1 (0+) 


— sin(7rap)- 


^ a p ,ot pfi- 1 


(A) 


— sin(7rap) 


00-/3,01/9+1 (^) 


sin(7rap)</i a p +1 , a p(A) + sin(nap)4>' a ~ ap+1 (0+) 


for A > 0. 


The two main results above, and in particular the techniques used to prove them, of¬ 
fer many new insights into the analysis of stable processes. Classically, the Wiener-Hopf 
factorisation of Levy processes provides the basis of many proofs for fluctuation identities, 
both exact and asymptotic. Historically there has been less exploration in this respect for 
the case of MAPs. But, nonetheless, the same importance of the role of the Wiener-Hopf 
factorisation applies, with many proofs following analogous lines of reasoning to the Levy 
case. See for example the Appendix in [[18]. When one now takes account of the degree 
of explicit detail that we offer here with regard to the Lamperti-stable MAP Wiener-Hopf 
factorisation, as well as the pathwise embedding of the fluctuations of this MAP into the 
fluctuations of stable process, one should expect to gain new results for the latter family of 
processes. Based on the computations derived in obtaining the matrix factorisations above, 
we offer some results in this respect at the end of this paper. Moreover, the robustness and 
applicability of the techniques we develop in proving the above two theorems also plays an 
important role in forthcoming work; see [28] and [27] . 


3 Non-symmetric Riesz—Bogdan—Zak transform 


A key component in proving Theorem [7] will be the use of the so-called Riesz-Bog dan-Zak 
transform which we now outline. 

Theorem 4 (Riesz-Bogdan-Zak transform). Suppose that X is a stable process as outlined 
in the introduction. Define 

rj(t) = inf{s > 0 : f \X u \~ 2a du > t}, t > 0. 

Jo 

Then, for all x 6 M\{0} ; (—1/X^(q) t > 0 under P x is equal in law to (X, P”-^), where 


dP° 


dP T 


Tt 


sin(7TQ!p) + sin(7ra/3) — (sin(7rap) — sin(7rap))sgn(X t ) A 
sin(7rap) + sin(7ra-p) — (sin(7rap) — sin(7ra-p))sgn(x) J 


X, 


x 


a—1 


lq<r{°>) (15) 


and Tt cr(X s : s < t), t > 0. Moreover, the process {X, P°), x G M\{0} is a self-similar 
Markov process with underlying MAP via the Lamperti-Kiu transform given by 


F°(z) = 


for Re(z) e (—a, 1). 


T(1 — z)T{a + z) 

T(1 — ap — z)T[ap + z) 

T(1 — z)T(a + z) 

r(ap)r(l — ap) 


T(1 — z)T{a + z) 

r(ap)r(l — ap) 

T(1 — z)T(a + z) 

T(1 — ap — z)T(ap + z) 


( 16 ) 
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In the case that X is a symmetric stable process (i.e. c + = c_, equivalently p — 1/2) 
the result is contained in the result of Bogdan and Zak [T2], who deal with isotropic stable 
processes in one or more dimensions. One may see the work of Bogdan and Zak, specifically 
the idea of spatial inversion through a sphere (here an interval), as building on original 
results of M. Reisz, who used this technique to analyse potentials, cf. [30] pp. 13-171], [3Tj 
as well as the discussion in Section 3 of [XT] , 

It is straightforward to deduce that (X, P°), x E R\{0} is a rssMp, inheriting the index 
of self-similar It iy a from (X, P x ), x E R\{0}. When a E (1, 2), |T6] have identified (X, P°) to 
be the law of a stable process conditioned to avoid the origin when issued from x G R\{0}. 
Their requirement that a G (1, 2) pertains to the fact that points are polar for a G (0, 1] and, 
accordingly, conditioning to avoid the origin makes no sense in the latter parameter regime. 
Nonetheless, the change of measure (fT7>D is still meaningful and gives preference to paths 
that approach the origin closely, penalising paths that wander far from the origin. In fact, 
we shall see in due course from its Lamperti-Kiu representation that (X, P°), x G R\{0}, is 
absorbed at the origin almost surely; see the forthcoming Remark [TJ In this sense, (A", P°), 
x G R\{0}, may be considered to be the stable process conditioned to be absorbed at the 
origin. When a — 1, one easily sees that P° = P x , x G R. 

In order to prove Theorem[4]we first need to briefly discuss the analogue of the exponential 
change of measure and Esscher transform for MAPs. Referring back to (fl3lh for each z E C 
such that Re(^) G (—l,a), there exists a leading real-valued eigenvalue of the matrix F(z), 
also called the Perron-Frobenius eigenvalue ; see j7[ §XI.2c] and [19] Proposition 2.12], If 
we denote this eigen value by x(z), then it turns out that it is larger than the real part of 
all its other eigenvalue. Furthermore, the corresponding right-eigenvector v(z ) has strictly 
positive entries, and can be normalised such that 7r • v(z) = 1, where we recall that tt is the 
stationary distribution of the underlying chain J. 

The leading eigenvalue x( z ) features in the following probabilistic result, which identifies 
a martingale (the analogue of the Wald martingale), a change of measure and the analogue 
of the Esscher transformation for exponents of Levy processes; cf. [7] Proposition XI.2.4, 
Theorem XIII.8.1]. (Note that the result is still true for general MAPs as introduced in 
Section IT21 1 

Proposition 3. Let Q t = a{(£(s), J(s )) : s < t}, t > 0, and 


M(t, 7) = e 7(€(t) ~ g(0)) ~ x(7)t l,J(t) [ 7 j , 

vj(o)(V 


t > 0 , 


(17) 


for some 7 such that x{l) defined. Then, M(-,y) is a unit-mean martingale with respect 
to (Gt)t> o- Moreover, under the change of measure 


dP^ 


dP, 


= M(t, 7 ), 


Qt 


t > 0 , 


the process (£, J) remains in the class of MAPs and, where defined, its characteristic exponent 
given by 

F-y(z) = A M^Fiz + 7) A„(7) - x( 7 )I, (18) 

where I is the identity matrix and A v (y) = diag(u( 7 )). (The latter matrix we understand 
to mean the diagonal matrix with entries of v(x) loaded on to its diagonal.) 
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Proof of Theorem [7} When we take F to be given by (fT3]h we can compute explicitly the 
quantity n as well a v(y) for a particular value of 7 that is of interest. We are interested 
in the case that gamma := a — 1. Note that 7 G (—l,a). A straightforward computation 
shows that, for Re(z) G (—l,a), 


det F(z) = 


r(a-z) 2 r(l + z) s 


7T- 


{sin(7r(«p — z)) sin(7r(ap — z)) — sin(7rap) sin(7r«p)} , 


which has a root at z — a — 1. In turn, this implies that xi a — 1) = 0. One also easily checks 
with the help of the reflection formula for gamma functions that 


v(a — 1) oc 


sin(7r«p) 

sin(7r«p) 


and, by considering -F(O) = Q, 


7T cx 


sin(7rap) 

sin(7rap) 


(19) 


We see that with 7 = a — 1, the change of measure (TT7T) corresponds precisely to (fl5l) when 
(£, J) is the MAP underlying the stable process. In particular, we can now say that the 
MAP associated to the process (AT, P°), x G M\{0}, formally named F°(z), is equal to 


F a -i(z) 


T(1 — z)T(ot + z) 

T(1 — ap — z)T(ap + z) 

T(1 — z)T(a + z) 

r(ap)r(l — ap) 


r(i — z)T{a. + z) 
r(ap)r(l — ap) 

T(1 — z)T{a + z) 

T (1 — ap — z)T{ap + z) 


for Re(z) G (—a, 1), where we have again used the reflection formula for the the gamma 
function to deal with the terms coming from A v (a — 1 ) in (fT 8 j) . 

Now let us turn our attention to the process (—1 /First note that, if X is an 
(a, p) stable process, then —X is a (a, p) stable process. Next, we show that (—1/A" r? (q ) t > 0 
is a rssMp with index a by analysing its Lampert-Kiu decomposition. 

To this end, note that, if (£*, J*) is the MAP that underlies A"* := —X, then its matrix 
exponent, say F*(z), is equal to (1T3]) with the roles of p and p interchanged. As X* is a 
rssMp, we have 


x t = ex P {£*0*0)) + i n(J*(v*(t)) + 1)} , t< r {0} , 


where 


rv*(t) 


i aC(s) ds = t. 


Noting that 


cvh) 


-2ar(^*(«)) du = ^(t) < r l°}, 


a straightforward differentiation of the last two integrals shows that, respectively, 

dif*(t) = e _ Qr(¥ ,*( t) ) and drift) _ or,®) ^ MO} 


dt 


d t 


= e 


pit) < 


12 















The chain rule now tells us that 


d(y?*Q y)(t) 
d t 


dy?*(s) 

ds 


s=r/(t) 


Mt) 

(it 


e a£*{<P*°v(t)) 


and hence, 


r<P*°V (*) 


3"“ S * (u) dM = t, 


T](t) < T®. 


The qualification that rj(t) < rf°i only matters when a 6 (1,2). In that case, the fact 
that P a ,(rl°l < oo) = 1 for all x G M implies that lim^oo^* = —oo almost surely. As a 
consequence, it follows that J 0 °° e~°^*^du = oo and hence lim^oo (p* o r](t) = oo. That is to 
say, we have lim^oo rj(t) = rl°l. Noting that for fceN, e~ mk = e 17rfc , it now follows that 


= ex P {-£>* ° V(t)) + i o ri(t)) + 1)} , t < r {0} 
v(t) 


is the representation of a rssMp whose underlying MAP has matrix exponent given by 
F*(—z), whenever it is well defined. Recalling the definition of F*(z), we see that the 
MAP that underlies (—l/X v ^))t>o via the Lamperti -Kin transform is identically equal in 
law to the MAP with matrix exponent F a _i(z). The proof is now complete. □ 


4 The ascending ladder MAP 


We shall derive the Matrix exponent k, by deriving each and every component of the matrices 
diag($i(A), <I> 2 (A)), A and K( A), for A > 0. In order to do this, we will make use of the 
Riesz-Bogdan-Zak transform from the previous section as well as some classical Markov 
additive renewal theory. In order to understand how the latter bears relevance, we need to 
briefly recall how the ascending ladder process ( H , J + ) emerges as a consequence of excursion 
theory and accordingly is a non-decreasing MAP. 

Let Y^ ] = (i V ((i)) - ((t), t > 0, where £(i) = sup s<t £(s), t > 0. Following ideas that 
are well known from the theory of Levy processes, it is straightforward to show that, as a 
pair, the process (Y^ x \ J ) is a strong Markov process. For convenience, write Y in place of 
F' ( °). We know by standard theory (c.f. Chapter IV of [ 18 ] ) there exists a local time of (V, J ) 
at the point (0,z), which we henceforth denote by {L[^ : t > 0}. Now consider the process 


L t :=J2 L t ] ’ t>0. 

iOE 


Note that this local time can be constructed uniquely up to a multiplicative constant, which 
will turn out to be of pertinence later. Since, almost surely, for each i ^ j in E, the points of 
increase of L W and are disjoint, it follows that (L~ 1 ,H + , J + ) := {(L^ 1 , H + (t), J + (t)) : 
t > 0} is a (possibly killed) Markov additive bivariate subordinator , where 

H+ (t) ■= £(4 _1 ) and J+ if) ■= J ( L i l )i if Li 1 < °°> 
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and H + (t) := oo and J + (t ) := f (a cemetery state) otherwise. Note, as a Markov ad¬ 
ditive subordinator, (L~ 1 , H + , J + ) is a Markov additive process with co-ordinatewise non¬ 
decreasing paths such that, in each of the states of J + , the process H evolves as a subordi¬ 
nator possibly killed at an independent and exponentially distributed time, whereupon it is 
sent to the cemetery state oo and J + is sent to f. Killing rates may depend on the state of 
J + . 

If we define 


e t = {e t (s) := + s) - : s < AL t 1 }, if A L t 1 > 0, 


and e t = d , some artificial isolated state, otherwise, then it turns out that the process 
{et : t > 0} is a (killed) Cox process. Henceforth, write n* for the intensity measure of 
this Cox process when the underlying modulating chain J + is in state i G E. As a Markov 
additive subordinator, the process (H + , J + ) has a matrix exponent given by 

Ec, ( J + (t) = j] = A > 0 , 

where k + (A) was given in (jUJ). Note in particular that, for i — 1 , 2 , $j(A) is the subordinator 
Bernstein exponent that describes the movement of H + when the modulating chain J + is in 
state i. Moreover, A is the intensity of J + and the matrix K{ A) = (K(X)) i j is such that, 
for i j in E, its (i, j)-th entry is the Laplace transform of the additional jump incurred by 
H when the modulating chain changes state from % to j. The diagonal elements of K( A) are 
set to unity. In general, we can write 

poo 

<Li(A) = rii(( = oc)+bi\+ (1 - e' Ax )nj(e^ G dx, J(() = i, ( < oo), A > 0, 

Jo 

where bj > 0 and ( = inf{s > 0 : e(s) > 0} for the canonical excursion e. For the case of the 
Lamperti-stable MAP, on account of the fact that the stable processes we consider in this 
paper do not creep, we can immediately set bj = 0 for i = 1,2. In the case that a G (0,1], 
points are polar and the underlying stable process explores arbitrarily large distances from the 
origin. When a G (1,2), the MAP in question represents the stable process until absorption 
at the origin, which occurs almost surely. We cannot rely on the range of the stable process 
until this time being unbounded and therefore 7ij(C = oo) > 0 , for i = 1 , 2 . 

Let us recall the following result, which is a special case of a general Markov additive 
renewal limit theorem given in the Appendix of [T 8 j. Such limit theorems are classical and 
can be found in many other contexts; see [29], [20J and 12, [3] to name but a few. First we 
fix some notation. For each a > 0, let 


T a = inf{t > 0 : H + (t ) > a}. 


Lemma 2. Suppose that the ladder height process ( H + , J + ) does not experience killing, that 
is to say, L 0 0 = oo. Then for x > 0 and i,j G {1, 2} ? 


lim P 0 ,j (H + (T a ) -aedx, J + {T a ) = j) 

a—>• oo 


l 

E 0 ,.(#+(!)) 


7 r i 7 l j( e (0 > x,J( 0 = j, c < oo) + 7 T fc A fc j(l - F fc + (x)) 


dx, 


( 20 ) 
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where k G {1,2} is such that k j and J, () ^ e Xx F^-{dx) = E[e~ xu ^}. /Is a more refined 
version of the above statement, we also have that 

lim P oAH + (T a ) - a G dx, J + {T a ) = j, J + (T a -) = j) 

a—^ oo 

= e 0 (H+(i)) njnj ( e ^ > x ’ = < °°) dx ’ ( 21 ) 

For all limits above, we interpret the right hand side as zero when Eo i7r (i/ + (l)) = oo. 

Recall that the local time L can be constructed up to an arbitrary multiplicative constant. 
If one follows how this constant permeates through to the definition of n*, i — 1, 2 and Ajj, 
i,j = 1, 2, we see that they are also defined up to the same multiplicative constant. For this 
reason, we shall, without loss of generality assume that this constant is chosen such that 
E 0)n (H+(l)) = 1. 

As will be explained in the forthcoming computations, the above Lemma provides the key 
to picking out the individual components that contribute to the matrices k(A) and k( A) by 
decomposing the left-hand side of (1201) and (12T| in terms of the behaviour of the associated 
stable process. Ultimately what we shall see is that all computations boil down to identities 
that come from the so-called two-sided exit problem for the stable process, which is originally 
due to [32] (see also Exercise 7.7 of [26]). 

To elaborate in a little more detail, let us introduce the stopping times for the stable 
process: for each a G M, 

rfi = inf{f > 0 : X t > a} and t~ = inf(t > 0 : X t < a}. 

Theorem 5. Suppose that X is a stable process as described in the introduction. Then, for 
9 > 0 and x G (0,1), 

P x (X r + — 1 g d0; 7i + < Tq ) 

= 8m ^ a P^ (X _ x) ap x a P9~ ap (e + 1 )~ ap (9 + 1 - x)- 1 ^. 
i r 

Equivalently, by scaling and translation, for 9 > 0 and x G (—1,1), 

P x (X Ti+ -lGd0;r 1 + <rr 1 ) 

= slll( ™ p ) (i _ x yp(\ + x) ap 9~ ap (0 + 2)~ ap (9 + 1 - x)~ 1 d9. 

7r 

Note, we have given two forms of the expression in the theorem above purely for convenience 
as they will both be used in the computations below. 


5 The ascending MAP with a £ (0, 1] 

When a G (0,1], we have that points are polar for the stable process X. In particular, as a 
rssMp, the origin is not accessible. Recalling that MAPs respect the same trichotomy as Levy 
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processes in terms of drifting and oscillating, we thus have that the Lamperti-stable MAP 
satisfies limsup^^ £(t) = oo. This means that the conditions of Lemma [2] are satisfied. 

In order to apply the aforesaid lemma, write X ( x ) to indicate the initial value of the 
stable process, i.e. Xq 1] = x E M\{0}. Thanks to self-similarity, and the Lamperti Kiu 
representation, we have, for example, that, on the event { J + (T a ) = 1}, i.e. {X T + Ar - > e“}, 


X {x } _ 

exp {H + (T a ) - a} = T ^ T -^ = d X {xe 


( 22 ) 


Hence, taking limits, we have, for i,j E {1, 2} and u > 0, 

lim P 0 ,i(H + (T a ) — a > u, J + {T a ) = 1) = lim P x (X T + > e u , r+ < tTJ. 

afoo x —>0 1 

Moreover, if we write X t = sup s<< X s and X_ t = inf, s < t X s , t > 0, then we also have 
7Ti77,i(e(C) > u,J( C) = 1,C < oo) 


= -T- limP x [X T+ > e u , X +_ > \X T+ _\^i < Tl, 
d u z-K) V r i T 1 1 r i 

= ~ f n x rf > e\X r f_ 6 dz,T* < rlj. 


(23) 


In order to progress our computations further and reach the goal of producing an identity 
for k(A), we shall first establish some intermediary results, starting with the following. 

Lemma 3. For A > 0, 

A sin(Trap) + sin(Trap) _ 

^ap-\-l,ap\^) • 

7r 

Proof. We start by noting that, for y E [0,1], 

P(W i+ > e\X r} _ < v.t* < r-) 

= P(X T j >e”,r+ <r-) 


= P 


z/(z+y) 


sin(7rap) 


71 


X + - 1 > 


y 


y 


z+y 


z + y 

Oip 


> r i < To 


y + z 


z + y 


Oip 


t~ ap {t + l)~ ap (t + 1 


z + y 


d t, 


where the penultimate probability follows from scaling and the final equality follows from 
Theorem [5j Hence, it follows that 

TTini(e(C) > «, 7(C) = 1, C < oo) 

= ~Lij-L P(A T > s y,E < 

= Sm(7rQ '^ [ e u ^-y ap z ap {e u - y)~ ap (e u + z)~ ap 

n Jo d y 

= 0l p Sm ^ K0t P) f e u z a ~ 1 {e u -z)- ap - 1 (e u + z )- ap dz 
71 Jo 


dz 


y=z 


(24) 
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Finally we can now take Laplace transforms and compute 


7Ti$i(A) = A 


= A ap 


e An 7Ti77,i(e(C) > u,J( C) = 1, C < oo) 
sin(7rap) 


1 /-oo 


= ap- 


7T 

sin(7rap) 



-(A+a)« z a-i( 1 + ze -»y°0(i _ Z e~ u )- {ap+1) du dz 


1 0 JO 

^OO 


7r 


(l-e- Aw ) 


(1 + e _u, ) Q ^(l - e - w ) a ^ +1 


du> 


sin(7rap) 


7r 


^ap+l,ap (A), 


( 25 ) 


where, in order to get the third equality, we have first substituted 6 = ze~ u in the second 
equality, then used Fubini’s Theorem and finally substituted 6 = e~ w . The result follows 
once we recall that 7Ti = sin(7rap)/(sin(7rctp) + sin(7rap)). □ 

We are now in a position to identify k( A) as presented in Theorem [2] when a G (0,1]. 

Proof of Theorem [H (i). In the spirit of Lemma [2], we also note that, again with the help of 
Theorem 0 we have, for u > 0, 

lim P o,i(H + (T a ) -a<u ; J + (T a ) = 1) 

a—>oo 

= P(W 1+ < e“;T+ < r“,) 

= Pj(W+ - 1 < |(e“ - 1 );t+ < Tq ) 
sin 


m(nap) fl\ h {e x) t _ ap ^ + + 1 / 2 )- 1 d f, 


7T 


(26) 


where we understand P 0i j be the law of the MAP with matrix exponent F issued from (0, i) 
Moreover, it follows that, for A > 0, 


©i(A) := lim P 0 Ae~ HH+{Ta) ~ a) -,J + (T a ) = 1) = 


sin(7 rap) 


—X u 


7T 


(1 - e"“)" p (l + e~ u ) a P 


rdn. 


Suppose now we define f(x) = e “ x (l — e x ) Qp (l + e x ) A straightforward compu¬ 
tation shows that 

ap ap 


f\x) = -f(x) 


+ 


(1 — e~ x ) (l + e-*) 

This will useful in the following computation, which also uses the conclusion of Lemmas [2] 
and El as well as integration by parts: 


TT2^2,lK(\)2p — @l(A) — 7Ti 

sin(7rap) 

71 


$i(A) 


(1 — e A:E ) , f°° (1 - e Ax ) , ap 

■fix) dx - / --- fix) j- -rdx 

Jo A Jy J (l- e~ x ) 


A 


sin(7rap) K ap ^p + i(A) 


7T 


A 
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for A > 0. 

Summarising the above computations, as well as the statement of Lemma El we have 
that, up to the multiplicative constant (sin (7r exp) + sin(7rcxp))/7r, 


$i(A) 


^ap-\-l,ap (A) 


and A 2 ,i-K'(A) 2 ,i = 


sin( 7 TQp) K a p } ap + i(A) 
sin (Trap) A 


A > 0, 


and hence, again up to the same multiplicative constant, 


sin(7rcxp) . 

A 2 .i = — 7 — 

smyKap) 


ap,ap +1 


( 0 +). 


By exchanging the roles of p and p, we similarly get expressions for <I> 2 (A), Ai^-K^Ajip and 
Ai ) 2 . Putting the pieces together into (EJ) we have the required form for k( A). □ 


6 The ascending MAP with a £ (1,2) 


When ex G (1, 2) the computations in the previous section break down as we must take 
account of the fact that the Lamperti-Kiu transform now only describes the stable process 
up to the first hitting of the origin. The way we will deal with this is to take advantage of a 
trick that emerges from the Esscher transform for MAPs. 

Revisiting Proposition [3] and the proof of the Riesz-Bogdan-Zak transform in Theorem 
[4) let us consider the MAP corresponding to (X, P°), x G M\{0}. Recall that its matrix 
exponent, F°(z), was given by (16). As well as being expressed via an Esscher transform of 
F(z) (see the proof of Theorem [jj, it also enjoys a Wiener Hopf factorisation so that 


— F°(id) = A~lk°(i Of 




{-id) 


= — A v (a — 1) 1 F(i6 + a — 1)A u (q; — 1), 


(27) 


for zGR, where n° is the stationary distribution associated to the underlying Markov chain 
of (X, P°), x G M\{0}, n° is the matrix exponent of the ascending ladder MAP of (X, P°), 
x G M\{0} and R° is that of its dual. It is easily verified that 


7T oc 


sin(7rap) 

sin(7rap) 


and hence, without loss of generality we may assume that n° = v(a — 1). Inserting the 
Wiener-Hopf factorisation into the expression on the right-hand side of (l27lh in a straight¬ 
forward fashion, one readily deduces that, for A > 0 


rc(A) = A 7 r oK°(A + a - l)A w l. 


(28) 


This can also be verified by directly performing the Esscher transform to the process ( H , J + ) 
through an application of Proposition [3] at the stopping time L t 1 ,t> 0. The reader should 
be careful to note that an additive shift of ex — 1 in the argument of F, the exponent of 
(£, J), corresponds to an additive shift of —(ex — 1) in the argument of k as the exponent of 
(H, J + ) is given by —k(—z). 
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Thanks to (l28]h it therefore follows that, to know k(A), it suffices to compute k°(A). 
This is favourable on account of the fact that the origin is polar for (X, P°), x G M\{0}, 
which implies that the MAP corresponding to F°(z ) does not drift to — oo and hence the 
conditions of Lemma [2] are met. As we shall soon see, in dealing with k°(A) via this lemma, 
we shall make effective use of the Riesz-Bogdan-Zak transform. Using obvious notation, let 
us write 

k°(A) = diag(<E>°(A), <f>°(A)) - A°o K°( A), A > 0. 

The analogue of Lemma El but now for k°, reads as follows. 


Lemma 4. We have 


^i(A) 


sin(7rap) + sin(7rap) 
sin(7rap) 


c(cr)0 a p_|_i 

, ap (A), 


A > 0. 


Proof. Referring to (l20]l and the discussion leading to (|23|) . we have, for u > 0, 


«(e(C) > u,J(C) = 1,C < oo) 

= —— lim / P° (x + > e u , X + G d z, \X + \ < z, rf < tZ 1 ) 
dux^oJ 0 x V T i r i “ T ! 1 V 

l /»1 

= —— limP° (x + > e u ,X + G dz, rjf < ru) 

d u J 0 o x \ T i r i “ 1 z ) 


d 

du 



, d m 

Inn —P 

dy 


o 

X 



< 


yxt 



dz. 


y=z 


Moreover, defining, for a < b, 


(29) 


= inf{t > 0 : X, 6 (a, 6)}, 

and writing P x , x G M, for the probabilities of —X (the negative of a stable process), we 
have, for 0 < x < y < 1 and u > 0, 


'imPx [X T + > X r+- < V, U < r_. 


du >o 


= _ i im P ° x r+ > e“, Ty < rZ z 

du x-X) x \ T v y 

= --^-limP_i /a .(JC T (-i/*,i/») e (—e — “, 0)) 

dn z—s-o 

= limPi /a .(X r (-i/, 1 i/ W ) G (0, e -u )) 

dn x-s-o 


P±oo 


f 2yze u — z + y A 2 yz 


{ 


y + z 


y + z 


( 30 ) 


where, momentarily, we will assume the limit, 


p ±00 (9) := lim P*(X T( -i,i) G d0)/dB, 9 G [-1,1], 

\z\—>oo 
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exists. In order to continue the process of identifying an explicit expression for (j30|) . we 
appeal a distributional identity found in [23] which gives the law of the position of a stable 
process when it first enters a finite interval. We can, in principle, take limits in the expression 
for IP^(A" r (-i,i) G d y), in particular, showing that p± is well defined. It turns out to be more 
convenient to root deeper into the proof of Theorem 1.1. of [25] and fish out an alternative 
expression. From equation (20) in [25) we have, for y G (—1,1) and a G (1,2), 


P+oJy) ■= lim P 2 (A r( -i,i) G dy)/dy 

z—> oo 

sin(7rap) 


7r 


-(1 + y)-°“>(l-y)-°“> 


x lim 

Z—>-oo 


(y + l)(z-l) c *(z + l)'+- 1 (z~y) 


-i 


+(1 - Dtp) T~ l / 3+1 t ap -\ 1 - tf~ a dt 


= (1 - ap)2 a ~ l Sm ^ ap \ l + y)~ ap (l - y)~ ap P t ap ~\ 1 - t) l ~ a dt 


7T 


r(ap)r(2 — a) 


^ ^ Y(2 — ap)Y(ap)Y(l — ap) 

— c(a)(l + y)~ ap (l — y)~ ap , 


(l + y)-'*(l-y)- c * 


(31) 


where 


c(a) = 2 


a—1 


T(2 — a) 


T(1 — ap)T(l — ap) 

By appealing to duality, one also easily verifies in a similar fashion that 


P-oo(y ) ■= lim P a (X T (-i,i) G dy)/dy = p^-y) = P +00 (y), 

z — y —oo 


where the function p^ is the same as p^ albeit the roles of p and p are interchanged. Hence 
we may accordingly refer to the function 

P±oc(y) ■■= P + oc(y) = P-oo(y) = c(a)(l + y)~ ap ( 1 - y)~ ap . (32) 


Returning to (1231) and (]30l) . we therefore have 

7r i ri i( e (C) >u,J( C) = 1, C < oo) = c(a) [ ^(ze~ u )e~ u dz, 

Jo 

where (3(9) = {p± oc (6 l )(6 l + 1) + p± oo (0)}/2c(a), which can easily be verified to satisfy 

r ( q \ = __ (ct — l)/2 

(1 - 9) a P +1 (l + 9)°p (1 - 6) a P(l + 9) a p’ 

for 9 G [—1,1]. We can now compute 


7r°$°(A) = c(a) A 


f* OO pi 


JO JO 

poo 


3 -(A+l )up( ze -u^ z d U 


= c(a) / (1 - e- Aw )/3(e- w )e~ w dw, A > 0, 
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as required. Note that the second equality condenses two changes of variable and an applica¬ 
tion of Fubini’s theorem, similar in spirit to earlier computations, into one step. The details 
are straightforward and left to the reader. □ 


Proof of Theorem El (%). Let Pg t be the law of the MAP whose matrix exponent is F°, issued 
from (0,z). Following previous reasoning, we compute, 

lim - a < u; J*(T a ) = 1) = lim K(X + < e«; T+< rl,) 

a—>-oo ’ x —>0 1 

= limP 1/x (X T (-i,i) G (e _u ,l)) 

>0 

= / P±oc(y)dy. 

Je~ u 

Hence, for A > 0, 

poo 

e°(A) := limP ° i (e- x ( H+ ™-^-J + (T a ) = l)= e-^ + d u p ±00 (e~ u )du, 

a-^oo ’ J Q 

and so, again appealing to obvious notation, from Lemma [2j 

,*i(A) 


(33) 


7r 2 °A°iT 0 (A) 2il = 0 ?(A)-tt°- 


A 


(1 - e _A “) f A/ 


A 


{p' ±00 (e~ u )e- u +p ±00 (e~ u )} e~ u du 


— c(a) 


(1 - e~ Xu ) 
A 


p(e- u )e~ u du 


(l - e _A “) 1 

A - j --5 {P'±oo(e"“)(e“" - 1) + P±oo(e-“)} e~“du 


= c(a )■ 


^ap,ap -\-1 


(A) 


A 


A > 0. 


This tells us that up to the multiplicative constant c(a)(sin(7rap)+sin(7rap))/sin(7rap) sin(7ra/3), 
for A > 0, 

X> ap,ap+ l(A) 


A 2 1 iT°(A) 2 ,i = sin(7rap)- 


A 


and A° ;1 = sin(7rap)^ p ^ +1 (0+). 


By exchanging the roles of p and p , we now have enough identities to complete fill out the 
entries of k(A). □ 


7 The ascending ladder MAP for the dual process 


Using (ITOll and (fTbl) a straightforward computation gives us 


F(z) 


r(« + z)r(i — z) 
F(ap + z)r(l — ap — z) 

r(« + z)r(i — z) 
r(ap)r(l — ap) 


r(« + ^)r(i — z) 
r(ap)r(l — ap) 

r(a + z)T{\ — z) 
Y(ap + z)r(l — ap — z) 


:i 
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which is well defined for Re(z) € (—a, 1). Glancing back, one quickly realises that F(z) 
takes the form of F a _±(z) (from the Proof of Theorem QJ, but with the roles of p and p 
interchanged; otherwise written 

F(z) = F„_,M \^ f = F°(z)\^ (34) 

for Re(z) G (—a, 1). This tells us that the dual of (£, J) is the MAP which underlies the 
rssMp (X, P°), x G M\{0} through the Lamperti-Kiu transform. 

We may now proceed to calculate £(A) by again taking advantage of the Esscher transform 
to remove the effects of killing. Indeed, referring to (I34p and noting that, without loss of 
generality, we may take 

V(a - 1)| pop = 7T, 

it is immediately clear that 

= A tt 1 «(* + 1 Ar, A > 0. 

Referring to (T28|) . this means, in particular that 

ft(A) = K°(A)| p ^, A>0. (35) 

That is to say £(A) is nothing more than k°(A) albeit with the roles of p and p interchanged. 

Remark 1. Whist (l35|) is helpful for computing £ in the case that a G (1,2), it also tells 
us something about n° in the case that a G (0,1), where k° has not already been evaluated. 
For the latter regime of a, appealing to Theorem [2] (ii), we note that £( 0 ) 1,1 > —£( 0 ) 1,2 and, 
similarly, £(0) 2 , 2 > —£( 0 ) 2 , 1 - This implies that the ascending ladder MAP associated to 
(X, P°), x G M\{0}, and hence (X, P°), x G M\{0}, is subject to killing. In turn, this means 
that, when a G (0,1), the MAP associated to F°(z ) drifts to — 00 and, accordingly, (X, P°), 
x G M.\{0}, is a real-valued self-similar Markov process which experiences absorption at the 
origin. 

In this sense, when a G (0,1), (X, P°), x G R\{0} may be reasonably named the stable 
process conditioned to be continuously absorbed at the origin. Indeed, in further work 
following ideas of [16] , we hope to give more mathematical substance to this remark. 

8 Cramer-type results for Lamperti-stable MAPs 

A classical computation that emerges from the Wiener-Hopf factorisation in the setting of a 
Levy process that drifts to — 00 , which also has a non-zero root of its characteristic exponent, 
is Cramer’s asymptotic estimate for the probability of first passage above a threshold, as 
well as the asymptotic conditional overshoot distribution conditional of first passage. See for 
example [9j. In the current setting we have noted that when a G (1, 2) and when a G (0,1) 
the ascending ladder processes of the Lamperti-stable MAP and and the ascending ladder 
processes of the dual Lamperti-stable MAP, respectively, undergo killing. Moreover, we have 
also noted the existence of roots to the leading eigenvalue of the associated matrix exponent. 
This means that we can expect to see Cramer-type results in each of these regimes. In this 
respect, we have two main theorems in this section. 
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Theorem 6. When a G (0,1], P o,i{T a < oo) = 1, for i = 1,2, where as, when a G (1, 2) 

T(2 -a) 


lim e ( “ 1)a P 0 ,i(T Q < oo) = 2 Q 1 — -——-- ——-rlp=i) + —- tv 

moo i (1 — ap)l (1 — ap) [sm(7rop) sm(7ra:p) 

Moreover, when a G (0, 2), for i — 1,2 and u > 0, 

lim P 0 ,i(H + T a ) — a G dw; J + (T a ) = j\T a < oo) 

a—>oo 

sin(7rap) , . . 

—-— —e~ au (l + e~ u )- ap {l - e~ u )~ ap du if j = 1 , 


L(i=2) 


= 


sm(7rap) e _ aU ( 1 + e _ u y ap ^ _ ^-ap^ ^ . = 2 


71 


Proof. The first claim follows by virtue of the fact that, as noted in Section [5j the ascending 
ladder MAP ( H , J + ) experiences no killing when a G (0,1]. Recall that 


Thanks to 


r ( i,i) i n f> 0 : X t G (—1,1)}. 
and Theorem [4l 


P0,1 (^a < OO) 

= P e -a(r+ A TZ 1 < r {0} ) 
= P)ie«(' r(_1 ’ 1 ^ < oo) 


= P°a(r ( 1,:l) < oo) 


(36) 


e —(a—l)a 

2 sin(7rap) 


E,, 


2sin(7rap)l (x (11)>0 ) + 2 sin(7rap)l ( x ( _ lil)<0 )) 


IQ;—1 


where, for each w G M\{0}, P° plays the role of P° with p and p interchanged (i.e. it plays 
the role of P^ for —X). Recalling the definition of the limiting distribution p± 00 given in 
), we thus have 

lim sin(7rap)e ( ' Q ” 1 ' )a Po i i(T a < oo) 

a—>oo 

= c(a) sin(Trap) [ y a ~ 1 ( 1 + y)~ ap { 1 - y)~ ap dy 


+c(a) sin(Trap) / y a (1 - y) ap (l + y) ap dy 
Jo 

, . sin( 7 rap) T(o;)r(l — ap) _ 

c(a)vr- —— --—2^i(ap, a, ®p + 1; -1) 

7r f (1 + ap) 

, . sm(nap) T(a)T(l — ap) H 

+c(a) 7 r-—- 1 _ I/1 -^- 2 -Pi(a/h a, ctp + 1 ; — 1 ) 


n T(1 + ap) 


= c(a) 7T, 


(37) 


where 2 -Pi (a, b, c; z) is the usual hypergeometric function and the final equality is a remarkable 
simplification which follows from one of the many identities for the aforesaid functions. See 
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for example the first formula at the functions.wolfrajn.com webpages [TJ. If, on the left- 
hand side of (157)) . we replace P 0 .i by P 0) 2 , the only thing that changes in the statement is 
that we must replace sin(map) by sin(map) on the left-hand side. This completes the proof 
of the first part of the theorem. 


For the next part, we split the proof into the cases that a G (0,1] and a G (1, 2). In the 
former case the result was already established in (1551) . For the latter case, appealing again 
to the Riesz-Bogdan-Zak transform, (1551) . (1571) and (1511) . we have, for u > 0, 


lim P 0 ,i (H + T a ) - a <u-, J + (T a ) = 1|T 0 < oo) 

a—> oo 

P e -a(X T+ < e“; t? < tZ 1 A t< 0 >) 


= lim 


a^-oo P c -a(r 1 + Ar_j < rf°l) 

P e 0 a (X r (- 1 , 1) G(e-“,l);r(- 1 ’ 1 )<oo) 


= lim 

a—>■ oo 


Pe“ ( r ^ -1,1 ) < °°) 
sin (map) 9 a ~ 1 p± 00 (9)d9 


sin(map) f 1 


m 


c(a)m 

9 a ~ 1 (l + 9)- a P(l-9)- ap d9, 


u > 0, 


(38) 


which is equivalent to the statement in the second part of the theorem when i — 1. It 
turns out that the asymptotic is unaffected when i = 2, however the details are left to 
the reader to verify. We also leave it as an exercise for the reader to check that when the 
event { J + {T a ) = 1} is replaced by { J + {T a ) = 2} on the left-hand side of (1551) . the resulting 
asymptotic is the same but with the roles of p and p interchanged. □ 


It is worth noting from the proof of this theorem that the methodology allows us access 
to new identities for stable processes with a G (1, 2). For example, the following polynomial 
asymptotic decay for the probability that the stable processes escapes (—1,1) on before 
hitting the origin. 

Corollary 1. Fora G (1,2), 


lim ^(a;)^ 1 "P x (r 1 + A r_ x < r®) 

x —>0 


i r (2 - tt ) 

T(1 — ap)T(l — ap), 


where 

. . sin(map) sin(map) 

s(x) := —^-^l(^>o) + —^-^l(*<o)- 

Another example of a new fluctuation result for stable processes is captured in the corol¬ 
lary immediately below. 

Corollary 2. For a G (1, 2), 9 > 0 and x G (0,1), 

F x (X r + - 1 G dO, r+ < tZ 1 A t {0} ) 

= Sm ^ ap ) M + x) al> (l - x) ap {2 + 9)~ ap 9~ ap { 1 + 9- x)- 1 
m 

-(a - 1 ) sm ( 7ra ^) ( 2 + 9)~ ap 9~ ap (l + 9)~ 1 x a ~ l [ (t - 1 ) ap ~ 1 (t + l) ap ~ l dt. 

TT J 1 
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Proof. According to Theorem [4j we thus have 


P x (X T + - 1 >e,T+< Tl 1 A r {0} ) 

= P x {—l/X T + G (-(1 + 0)-\ 0), r+ < rA A t<°>) 

= P°_ 1/3; (X t( - M) € (-(1 + 0)~\ 0), < oo) 

= x a - l t l/x (\X T ^ (“- 1 ; X r( - M ) G (0, (1 + 0)" 1 )), 

for 6 > 0. It follows that 

W x {X T + - 1 G d0, r+ < TZ l a r {0} ) = x*- 1 &/* {iTe) d0 ’ 

where, for x > 1 and y G [—1,1], p x (y) = P x (X r (-i,i) G dy)/dy. The latter can be found in 
Theorem 1.1 of j25| and is given by 

P,(A' r , e dy)/dy = S,n(7rap) 0 + 1 r(x - 1)^(1 + 1 - v)-°*(x - y)- 1 

7r 

- (a - 1 ) sm ( 7ra P) ( 1 + y )-“P(i _ f\ t - l) ap -\t + l)^- 1 dt. 

^ J i 

The result now follows by straightforward algebra. □ 


Let us now turn our attention to the dual of the Lamperti-stable MAP when a G (0,1). 
We denote its law by P,,.^, for x G M and i — 1,2. Recall from Section [7] that this MAP 
corresponds to the rssMp (X , P°), x G M\{0}. The analogue of Theorem [6]takes the following 
form. 


Theorem 7. If a G (0,1), then, for i = 1,2, 
lim e (1 -“ )a P 0li (T a < oo) = 


T(ap)T(ap)T(2 — a) (sin(7rap) 


7T 71 

-l(»=i) H—:—7-—1 


sin(7rap) 


(*= 2 ) 


If a G [1,2), then P o,i(T a < oo) = 1, fori = 1,2. Moreover, for a G (0,2), i = 1,2 and 
u > 0, 

lim P oi(H + T a ) - aGdii; J + (T a ) = j|T a < oo) 

a—>oo 

T(2 — a) 


')OL— 1 


= 


')(%.— 1 


r(l — ap)T(l — ap) 
r(2 — a) 


T(1 — ap)T(l — ap) 

Proof. Appealing to (j22|) . we start by noting that 


e"“(l + e _ “) _ “P(l - e~ u )~ ap du if j = 1 
e"“(l + e _ “) _ “P(l - e~ u )~ ap du if j = 2. 


Po,i(T a < oo) = Pg—a( t^ - A t_ x < oo), 
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and hence, by again making use of Theorem [5] and Theorem [2 we also have that, when 
a G (0,1), 


lim sin(7r«p)e^ a ^°Po,i(T a < oo) 

a—>oo 

= lim E e -a ( ( sin(7rap)l ( x _ >o) + sin(7rap)l( X + _ <o) ) \X +^ ~ 

a —>oo \ V T Ar -i t i At -i ' 1 _1 

3 (a-l)«^« _ + l)-“Pd M 

XD 

e ( Q -i)“( e u - l) _ap (e“ + l)“ a/5 d u 


IQ —1 


sin(7rap) sin(7rap) 




n 

Jo 

sin( 

nap) sini 

{nap) 



n 


sin( 

[nap) sin(7rap) 

7r(l 

— 

ap){ 1 - 

ap) 

sin( 

nap) sin 

{nap) 

7T 

(1 

- ap){l 

- ap) 

sin( 

nap) sin(7rap) 

7r(l 

— 

ap){ 1 - 

ap) ' 



2 1 -“tt 



T{ap)T{ap)T{2 — a) 


(1 - ap) 2 J r i(l, ap, 2 - ap ; -1) 


T(2 — a) 


(39) 


where the penultimate equality is again remarkably due to a very particular identity for 
hypergeometric functions; see the second formula in the functions.wolfram.com webpage 
[1 j. If we repeat the computation with P 0) i replaced by P 0i2 , then the only other thing that 
changes in (l30l) is that sin(7rap) is replaced by sin(7rap) on the left-hand side. This completes 
the proof of the first part of the theorem. 

When a G [1,2), the ascending ladder height MAP of the dual is not killed (see the 
discussion in in Section [7]) and hence P 0 ) i(T a < oo) = 1, for i — 1,2. 


For the next part set a G (0,1]. Starting as we did in the proof of Theorem [2] (i), we note 
from Lemma [2j (39) and Theorem [5] that 

lim P 0 ,i {H + (T a ) — a < u; J + (T a ) = 1| T a < oo) 

a—> oo 

= lim P°- a (A" r + < e u ; rf < rZ 1 |r 1 f ' A tZ 1 < oo) 


= lim 

a—>oo 


P° e _ a (X T +- 1 < e u - 1 ,r 1 + <rr 1 ) 


P°-a(r 1 + A T_! < oo) 


r(ap)r(2-a) A x 

- T(1 — ap) nX rt -V- 

= 2 «-i r(2 - a) p '- 1 

T(1 - ap)T(l - ap) J 0 


l<e u -l ,r+<r: 1 ) 

(t9 + l) a - 2 6- a P(6 + 2)~ ap d 6. 


This is equivalent to the second statement of the theorem for i = 1. The computation when 
i — 2 can be performed similarly. Replacing the event { J + (T a ) = 1} by {J + (T a ) = 2} in the 
probability above, affects the final equality only exchanging the roles of p and p. The details 
are left to the reader. 
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Filially, when a G (1, 2), the desired asymptotic can already be found in ([33]) . as soon as 
one notes that P 0 ,i agrees with Pq - when the roles of p and p are exchanged. □ 

Similarly to before, one can proceed to extract further identities for the Doob h-transformed 
process (X, P°), x G M\{0}, however, we leave this for the reader to amuse themselves with. 
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